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Mechanism of Stepped Leaders in a Simple Discharge Model 
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We construct a one-dimensional model for the stepped leader in the filamental discharge by 
simplifying an electric-circuit model of discharge. We find that the leader of the discharge moves 
stepwise by direct numerical simulations, and then we try to understand the mechanism of the 
stepwise motion by reducing the spatially extended system to the dynamics of the tip position of 
the discharge. 
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I. INTRODUCTION 



The discharge occurs when the voltage between two electrodes is beyond a critical value. There are various forms 

of discharge, such as corona discharge, spark discharge, glow discharge, and arc discharge, depending on various 

,^ \ conditions, such as pressure, temperature and the shape of electrodes. For gas discharge, the pressure p and the gap 

Q-c length I between the electrodes are important. When p x Z is large, the spark discharge takes a form of a filament. 

^ ■ Meek proposed the streamer theory for the filamental discharge. [l| The tip region of a growing filamental discharge 

is sometimes called a leader. It is known that the leader grows stepwise in the lightning discharge. The gap distance 

between the electrodes is very long in the lightning discharge. Ql The leader goes down from a thunder cloud to the 

ground relatively slowly. The velocity is around 10^ m/s. The extending process of the leader is invisible to the 

eyes, but it can be observed with a high-speed camera. The leader moves in steps of about 30 m with a pause of 



^ \ about 40 ms between steps. This type of leader is called the stepped leader. When the leader reaches the ground, 

^ • a strong flash called a return stroke appears instantly, which we observe as lightning. However, the origin of the 

£^ ' stepped motion is not completely understood. There were several qualitative theories for the stepped leaders. Bruce 

pointed out the importance of the transition from a weak glow discharge to a strong arc discharge. [3] Kumar and 

Nagabhushana proposed a simulation model of stepped leaders based on a complex electric breakdown process. [^ In 

a previous study, we proposed a simple deterministic electric circuit model composed of resistors and capacitors, and 

performed a numerical simulation on triangular lattices to iirvestigate complex patterns in discharge processes. fSj In 

K^ \ the model, a two-step function for the conductance is assumed, which expresses a transition from a weak discharge to 

, ' a strong discharge. Branched patterns of the discharge similar to the Lichtenberg figure were found in the numerical 

,_^ \ simulation of the model. We found a stepwise motion of the leaders in a certain parameter range and showed that the 

fSl branching of the pattern is not directly related to the stepped motion. In this paper, we simplify the electric circuit 

^^ model to a one-dimensional model and try to understand a mechanism of the stepped motion qualitatively. 

o 

II. ONE-DIMENSIONAL MODEL AND NUMERICAL SIMULATION 



We consider a one-dimensional electric circuit model in this paper to simplify the argument. The gap length 
between the electrodes is assumed to be L, and the voltage Vq is applied between x = Q and x — L. In the numerical 
simulation, the space of size L is discretized with the interval Ax. The interval Ax = 0.1 is used in most numerical 
2 . simulations. The electric potential at the ith site is denoted by Vi at i = x/Ax. A resistor is set between the ith site 
■ ■ ' and the {i + l)th site, and the conductance is expressed as ai/Ax. A capacitor is set between the ith site and the 
earth, and the capacitance is assumed to be C • Ax. Here, we assume that the conductance is inversely proportional 
to the interval Aa; and the capacitance is proportional to Ax. The time evolution of Vi is expressed as 

C^ = {a,_i(F,_i - V,) - adV, - l/,+i)}/(Ax)2. (1) 

By the continuum approximation, eq. (1) is reduced to the partial differential equation 

dV d ( dV\ 

where V{x) and (j{x) denote Vi and ctj at i = x/ Ax, respectively. The boundary conditions for V{x) are expressed 
as V{x) = Vo at a; = and V{x) = at x = L. Similarly to the previous paper, we assume a two-step function for 
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FIG. 1: (a) Time evolution of the tip position xi{t) of the discharge for Vo = 1,C — 0.03, Aa; — 0.1, a — 0.02, ai — 50, (T2 = 
2000, r = 0.01, Scio = 0.1,£c20 = 10, and L = 90. The dashed line denotes xi{t) = vit + xo with «i = 280. (b) Snapshot 
profile of a-{z) where z = x ~ x\. The dashed curve denotes the theoretical curve obtained using eq. (5). (c) Snapshot profile 
of V{z) in the range of —5 < 2 < 0. The dashed line denotes V = —0.1 • z. (d) Snapshot of Viz) in the range of —22 < z <\. 
The solid curve denotes the numerical result and the dashed curve denotes eq. (7). 



the conductance, although we have observed a similar behavior even if the two-step function is slightly modified to a 
continuous function. Thus, the time evolution of Oi is expressed as 



~dt 

~dt 
dfJi 

~dt 



-a, for E, < Eel, 

a I - a, for Ed < -B, < Ec2, 

(72 - a, for Ec2 > Ei, 



(3) 



where Ei = \Vi+i — Vi\/Ax is the local voltage between the ith and {i + l)th sites, and Ed and Ec2 are threshold 
values for the weak and strong discharges, t is the relaxation time, and cri^2 are stationary values of the conductance. 
We assumed that (72 is much larger than cti because the ionization proceeds rapidly at the transition from the 
weak discharge to the strong discharge. We further assume that the threshold Ed and Ec2 decrease with a as 
Ed = Edo ~ eta and Ec2 = £-020 ~ cua. If the local voltage is below the first threshold Ed in the entire region, 
the conductance decays to 0, which implies the insulator. When the local voltage is increased and goes beyond the 
first threshold, the discharge occurs and the conductance becomes nonzero. Then, the first threshold value slightly 
decreases as Ed = Edo — cua. This effect induces a hysteresis in the conductance. That is, if the discharge occurs 
once, the discharged state is maintained even if the local voltage is slightly decreased. A similar hysteresis is assumed 
to occur owing to the term Ec2 = £-020 — ota at the transition from the weak discharge to the strong discharge. This 
type of hysteresis is often observed even in experiments of discharge phenomena. 

If Ec2 is sufficiently large, the strong discharge does not occur, because Ei does not reach the second threshold. In 
this case, the leader or the tip position of the weakly discharged state moves smoothly. We show a numerical result 
in Fig. 1. The parameter values are Vq = \,L^ 90, Ax = 0.1, C = 0.03, t = 0.01, a = 0.02, cti = 50, ^2 = 2000,^:^1 = 
0.1, and Ec2 = 10, and the initial conditions are set to be (Xi = and Vi = 0. Figure 1(a) shows the time evolution 
of the tip position of the weak discharge X\ = iiAx. The tip position xi increases monotonically, and the velocity 
vi = dxi/dt around t = 0.08 is evaluated at vi ~ 280. However, the propagation velocity slightly decreases in time. 
This is because the distance between the electrode at j = and the tip position ii increases in time. Figures 1(b) 
and 1(c) show profiles of the conductance a{z) and the electric potential V{z) at i = 0.08. Here, z ~ x — xi denotes 
the distance from the tip position xi. For z > 0, a{z) = and V{z) = 0, which implies that z > {i > ii) is an 
insulator region. The conductance (t{z) increases gradually to ui as \z\ {z < 0) is distant from z = 0. The dashed 
line in Fig. 1(c) denotes V{z) — —Edoz, which implies that the local voltage E at the tip of the discharge is equal to 
Edo- Although the local voltage E = \dV/dz\ is smaller than Edo in almost the entire region of z < except for the 
region z ^ 0, the discharged state is maintained because of the hysteresis effect, that is, E is smaller than Edo but 
larger than Ed = Edo — aa{z). 

If we assume that the tip position Xi moves at a constant velocity vi, steadily moving solutions a{z) = (t{x — vit) 
and V{z) = V{x — vit) can be obtained from eqs. (2) and (3). Equation (3) leads to 
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FIG. 2: (a) Time evolutions of the tip positions x\(i) and X2(t) of the weak and strong discharges for Vq = 1, Ax = 0.1, C = 
0.03, a = 0.02, (Ti = 50, 0-2 = 2000, r = 0.01, Bcio = 0.1, Ec20 = 1, and L = 90. (b) Time evolutions of the tip positions X2{t) at 
Ec20 = 0.8, 0.9, 1, 1.1, and 1.2 for Vo — 1. (c) Time evolutions of the tip positions X2{t) at Vo — 0.3, 0.7, and 1 for i5c20 ~ 1- 
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FIG. 3: (a) Snapshot profiles of (^{z) where z — x — xi. The dashed curve denotes the theoretical curve obtained using eq. (5). 
(b) Snapshot profile of V{z) at t = 0.1 where z = x — xi in a. stepped stage. The dashed curve denotes a piecewise linear 
function of x. (c) Snapshot profile of V{z) a.t t — 0.09385 where z = x — xi in a moving stage. The dashed curve denotes a 
linked curve of a piecewise linear function (10) and the function expressed by eq. (11). (d) Snapshot profiles of (j{z) at t — 0.1 
(solid curve) and t — 0.09385 (dashed curve). 



for z < 0. The solution to eq. (4) is expressed as 

ct{z) = cri[l - exp{z/(t;ir)}]. 



(5) 



for z < and cr(z) = for z > 0. The velocity vi is evaluated at 280 near t = 0.08 from Fig. 1(a). The dashed curve 
in Fig. 1(b) denotes eq. (5) using vi — 280, r — 0.01, and cti = 50. Good agreement is seen between the numerical 
and theoretical curves. On the other hand, eq. (2) is reduced to 



dV 
-v^CV^aiz)^. 
oz 

The substitution of eq. (5) into eq. (6) yields the solution 

V{z) = V{zq) exp[(-«2Cr/fTi){(z - zo)/(z;ir) + ln(l - e'^''^-'^^) 



ln(l 



z/{vit] 



)}], 



(6) 



(7) 



where zq is a certain position where V{z) — V{zq) is satisfied. Figure 1(d) shows a comparison of V{z) (solid curve) 
obtained by a direct numerical simulation and V{z) (dashed curve) expressed by eq. (7), where V{zq) = Vq = 1 at 
zq = —20.1 (or x = 0) is used. The deviation of the two curves is considered to originate from the nonstationarity 
of this system, that is, the boundary condition V{zq) = Vq at the left electrode 29 — —vit is a moving boundary 
condition in this moving frame. 

We have performed several numerical simulations by changing the second threshold Ec2o &t a fixed value of Vq = 1 . 
The other parameter values and the initial conditions are the same as those previously used. We have found that 
the strong discharge does not appear at £"£20 > 1-8. The strong discharge occurs for Ec2o < 1-8, and the tip position 
X2 of the strong discharge moves stepwise for 0.85 < Ec2o < 1-17. X2 moves smoothly for Ec2o < 0.85. Figure 2(a) 
shows time evolutions of the tip positions xi and X2 respectively for the weak and strong discharges at Ec2o = 1 and 
Vq = 1. It is clearly seen that X2{t) moves stepwise, that is, X2{t) repeats the forward motion and the stop. On the 



other hand, xi{t) moves smoothly. Figure 2(b) shows time evolutions of X2{t) for £^£20 = 0.8, 0.9, 1, 1.1, and 1.2. The 
average velocity of X2{t) decreases with i?c20- The period of the stepped motion is the shortest near -Ec20 = 1- 

Several numerical simulations were also performed by changing Vq at a fixed value of Ec2q ~ 1- Figure 2(c) shows 
the time evolution of X2{t) for Vq =- 0.3,0.7, and 1. Stepwise motions are clearly observed. The average velocity of 
the tip position slightly increases with Vq; however, the period of the stepped motion changes rather strongly with Vq. 

The stepped motion of the leader has been studied in more detail for -Ec20 = 1 and Vb = 1. The tip position X2 
of the strong discharge exhibits a stepped motion, but the tip position x\ of the weak discharge moves smoothly, as 
shown in Fig. 2(a). The velocity of the first tip x\ of the weak discharge is evaluated as wi = 364 at the parameter 
values. Figure 3(a) displays several snapshots of (t(z) at different times, where z is the distance z — x — x\ from the 
first tip x-y. The dashed curve is u{z) = a\\V — exp{z/(wiT)}]. In the region of the weak discharge, the stationary 
solution of cr(z) by eq. (5) is rather good approximation. When the stronger discharge sets in, (t(z) increases rapidly. 
Almost discontinuously jumped lines in Fig. 3(a) correspond to the transition to the strong discharge. 

Figure 3(b) displays a snapshot pattern of V{z) at i = 0.1 when the leader is in a stepped stage. The position of X2 
remains constant with time in the stepped stage. The conductance a{z) at t = 0.1 is shown in Fig. 3(d) with a solid 
curve. The conductance is very large, i.e., cr ~ (72 ^ "'2 for z < Z2 = X2 — x\^ and rather small, i.e., ct ~ crj. ~ '''i for 
Z2 < z <^. If a'^ and a'2 are assumed to be certain constant values, V[x) can be roughly approximated at a piecewise 
linear function: 

V{x) - Vb + Ili_l^, for < x < X2, (8) 

X2 
V2{xx-x) c ^ ^ fr.-. 

— -, tor 2:2 < x < xi, (9) 



Xx - X2 



where 



Y2 = V{X2) ^° 



1 -^o'-yX2l{a'2{xx -X2)Y 

If o'y = 40 and a'2 = 1200 are used, V2 is evaluated at V2 = 0.811 for X2 = 37.7 and xi — 44 corresponding to the 
snapshot profile shown in Fig. 3(b). The dashed curve in Fig. 3(b) denotes this piecewise linear approximation of 
V{z). 

Figure 3(c) displays a snapshot pattern of V{z) at t — 0.09385 when the leader is in a moving stage. The second 
tip X2 moves approximately with the velocity V2 ^ 1640 ^ 4.5wi. The shoulderlike structure is characteristic of 
the profile of V{z) in contrast to the profile in Fig. 3(b), in which the local voltage E = \dV/dz\ is rather large 
near z — Z2. The conductance <j{z) is shown in Fig. 3(d) with a dashed curve. The conductance cr{z) and the 
derivative of V{z) have a discontinuity at 2:20 ~ —7.8. Here, the discontinuity point X20 = xi + Z20 is the position 
where the tip X2 of the strong discharge remained in the previous stepped stage. The conductance in the region of 
X2 < a; < xi is approximated at cri[l — exp{(x — xi)/{vit)}] and the conductance for X20 < x < X2 is approximated 
at CTi[l — exp{(x2 — xi)/{vit)}] + {02 — o'i)[l — cxp{(x — X2)/(w2t)}]. If the conductance is roughly approximated at 
certain constant values as a{x) = CT2 for < x < X20 and a(x) = a'^ for X20 < x < X2, Vix) is roughly approximated 
at a piecewise linear function: 

^(^) = y^+^Y2^^1^ for 0<x<X2o, 
X20 

T. , (V20 - V"2)(X20 -x) 

= V^o H forx2o<x<X2. (10) 

X2 - X20 

On the other hand, we use eq. (7) for V(x) in the region of X2 < x < xi: 

V{x) ^ V{X2) CXp[(-w2CT/ai){(x - X2)/(wit) + ln(l - e(^2-:ri)/(«ir)) _ l^(^ _ e(=r-=ri)/(fir))}] 

for X2 < X < x\. (11) 

The parameters V20 = V^(x2o) and V2 = 'V{x2) are unknown. They are determined by the boundary conditions of 
V{x) at X = X20 and X2 expressed by 

, Vb ~ V^20 , T/20 - V2 , V^20 - V^2 , ^ dV{x2) ,.^. 

0-2 =(^1 > CTl = -Cr X2^^ , 12 
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FIG. 4: (a) Time evolutions oi aa{x2),E{x2) and E + aa obtained using eqs. (14)-(17). (b) Comparison of the time evolutions 
oi a<j{x2) and E{x2) obtained using eqs. (14)-(17) (solid curves) and the direct numerical simulation (dashed curves), (c) Time 
evolutions of the tip position X2{t) of the strong discharge. The dashed curve denotes the numerical result and the solid one 
denotes the theoretical approximation. 

which are due to the continuity of the current. The parameters T^o and V2 are explicitly expressed as 



T^20 

V2 



o-^{cri + viC{x2 - a;2o)}Vb 



O"lV20 

a[ +viC{x2 - 2:20) 



2^20 )} +(t[viC'x2o' 



(13) 



The dashed curve in Fig. 3(c) shows a linked curve of the piecewise linear function (10) and the function expressed 
by eq. (11) using (T2 = 1800 and a[ = 100, although the horizontal axis is shifted to z — x ~ xi. We have not 
yet succeeded in obtaining the solution satisfying the correct moving boundary conditions with different velocities 
dx2a/dt — 0, dx2/dt — V2, and dxi/dt — vi. However, the above approximate solution reproduces a shoulder structure 
and a sharp derivative near x = X2 fairly well. 

III. SIMPLE MODEL FOR STEPPED LEADER 

On the basis of these numerical observations and the rough approximation for the profile of the electric potential 
shown in the previous section, we propose a very simple model for the time evolution of the position X2{t). In the 
stepped stage, the first tip xi moves with the velocity vi as a;i = xq + vit, where xq is a certain initial position, and 
the second tip X2 remains at a certain position: X2{t) = X20. In this stepped stage, the conductance a{x2) at a; = X2 
increases as 



cr{x2) = (Ti [1 - exp{-(xo + vit - X2o)/{vit)}], 



(14) 



because the difference xi — X2 — xq + vit — 2:20 increases with time. On the other hand, the derivative of the electric 
potential at x = X2 is evaluated using eq. (9) as 



E- 



a'Vo 



a[x2o + cr2(a;i -2:20)' 



(15) 



where a'l = 40 and a'2 = 1200 are used in the following numerical simulation. In the stepped stage, E decreases 
monotonically with time, because the difference Xi — X2 = xq + Vit ~ a;2o increases with time. The time evolutions 
of cr(2:2), E{x2), and E{x2) + acr(2;2) are shown in Fig. 4(a). The sum E + aa{x2) increases monotonically with time 
and reaches the second threshold Ec2o — 1 from below. Namely, the local voltage E{x2) goes beyond the threshold 
Ec20 — aa{x2)- Then, the gate to the strong discharge opens and the second tip 2:2 starts to move with the velocity 

V2- 

In the moving stage, the time evolution of X2{t) is expressed with 2:2 — 2:20 + V2{t — tn), and the conductance ct{x2) 
is approximated using eq. (5) as 



(t(2;2) = cri[l -exp{-(2;o + vit - X20 - V2{t - t„))/(i'ir)}]. 



(16) 



Here, i„ is the time when a transition to the moving stage occurs. The conductance a{x2) decreases monotonically 
with time in the moving stage because the distance xi — X2 — xg — 0:20 + V2tn — (f2 — Vi)t decreases with time. The 
derivative E — dV/dx of the electric potential at a; = X2 is evaluated using eq. (11) as 

{<7[ + viC{x2 - a;2o)}o'(a;2)-yiC" 

where the parameter values (t[ = 100 and (72 = 1800 are used. There is a discontinuity in E^s of eqs. (15) and (17) at 
i = i„ in this simplified model. The summation of E{x2) +017(3:2) decreases monotonically with time and reaches the 
threshold £"020 = 1 from above at i = i^ > tn- Then, the moving stage changes into the stepped stage because the local 
voltage E is below the second threshold i?c20 — ota, and X2 stops. The stepped stage continues again until t — i„+i. 
The repetition of the stepped and moving stages reproduces the behavior of the stepped leader, as shown in Fig. 4(c). 
Figure 4(b) shows time evolutions of a{x2) and E(x2) in the direct numerical simulation (dashed curves) and the 
theoretical approximation (solid curves). The time evolutions of a{x2) and E{x2) in the direct numerical simulation 
are qualitatively similar to the theoretical approximation. Figure 4(c) shows a comparison of the time evolutions of 
X2{t) in the direct numerical simulation (dashed curve) and the theoretical approximation (solid curve). The period 
of the stepped motion in the direct numerical simulation is the same order but slightly larger than the theoretical 
model. We think that the main reason for the deviation is the rough approximation obtained using eqs. (15) and (17) 
for the derivative of the electric potential dX x — X2- 

IV. SUMMARY 

We have proposed a one-dimensional model for stepped leaders and a simplified dynamical model for the position X2 
of the stepped leader to understand the stepwise motion. It was shown that the two-step function of the conductance 
and the shift of the threshold by the form Ec2 = E(.2q — ota is essentially important for the stepwise motion in our 
model. Our model might be very simple for the realistic discharge process, but we think that such a simplified model 
is useful for understandingthe unique motion of the stepped leader. Stepwise growth was observed in other systems, 
such as bacteria colonies. i(l| Some similar mechanisms might work also in these systems. 
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